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ABSTRACT 

We present the results of weak gravitational lensing statistics in four different 
cosmological ./V-body simulations. The data has been generated using an algorithm 
for the three-dimensional shear, which makes use of a variable softening facility for 
the iV-body particle masses, and enables a physical interpretation for the large-scale 
structure to be made. Working in three-dimensions also allows the correct use of the 
appropriate angular diameter distances. 

Our results are presented on the basis of the filled beam approximation in view 
of the variable particle softening scheme in our algorithm. The importance of the 
smoothness of matter in the universe for the weak lensing results is discussed in some 
detail. 

The low density cosmology with a cosmological constant appears to give the broad- 
est distributions for all the statistics computed for sources at high redshifts. In par- 
ticular, the range in magnification values for this cosmology has implications for the 
determination of the cosmological parameters from high-redshift Type la Supernova?. 
The possibility of determining the density parameter from the non-Gaussianity in the 
probability distribution for the convergence is discussed. 

Key words: Galaxies: clustering — Cosmology: miscellaneous — Cosmology: grav- 
itational lensing — Methods: numerical — Large-scale structure of Universe 



1 INTRODUCTION 
1.1 Outline 

We present the results of a study of the weak gravitational 
lensing of light in four different cosmological models, using 
the algorithm for the three-dimensional shear developed by 
Couchman, Barber and Thomas (1999). Since weak lensing 
effects depend on the angular diameter distances for the 
lenses and sources, and also the specific distribution and 
evolution of matter, the results are sensitive to the particular 
cosmological model. 

In strong lensing studies, frequent use is made of the 
'thin-screen approximation,' in which the mass distribution 
of the lens is projected along the line of sight and replaced 
by a mass sheet with the appropriate surface density profile. 
Deflections of the light from the source are then considered 
to take place only within the plane of the mass sheet, making 
computations for the light deflections simple. The simplicity 
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of the thin-screen approximation has also lead to its frequent 
use in weak gravitational lensing studies, where each of the 
output volumes from cosmological iV-body simulations is 
treated as a planar projection of the particle distribution 
within it. 

However, a number of problems can arise with two- 
dimensional approaches, especially in weak lensing studies 
in which the large-scale distribution of matter extending to 
high redshifts is responsible for the lensing. Couchman, Bar- 
ber and Thomas (1999) considered some of the shortcom- 
ings and were motivated to develop an algorithm to com- 
pute the six independent components in three dimensions 
of the second derivative of the gravitational potential (the 
three-dimensional shear). In the method used in the present 
work, for which the effects of lensing along lines of sight are 
required, the two-dimensional 'effective lensing potentials' 
(see Section 2) are obtained by integrating the computed 
three-dimensional shear components along the lines of sight. 
These effective lensing potentials are used to construct the 
Jacobian matrices which are recursively generated along the 
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lines of sight, and from the final Jacobian matrices the mag- 
nifications and two-dimensional shear are determined. 

A brief outline of this paper is as follows. 

In Section 1.2 we summarise the previous weak grav- 
itational lensing methods by other authors. 

In Section 2, the essential equations for gravitational 
lensing are explained, and the multiple lens-plane theory is 
described, culminating in expressions for the magnification, 
convergence, two-dimensional shear, and ellipticity, which 
are the weak lensing outputs from the cosmological simula- 
tions. 

Section 3.1 summarises the three-dimensional shear 
algorithm, and particularly the choice of the variable soft- 
ening scale for the particles. This feature limits the effects of 
isolated particles (not representative of the large-scale struc- 
ture), and builds in a physical significance both to the choice 
of the softening scale, and to the underlying form of the dark 
matter. We explain also how the softening carefully limits 
the incidence of strong lensing events in all the cosmologies. 
In Section 3.2 we describe the four cosmological iV-body 
simulations used in this work. In Section 3.3 we state how 
the simulation boxes are oriented to minimise correlations 
in the large-scale structure between adjacent time-outputs, 
and we describe the establishment of lines of sight through 
the simulations, and the locations of the evaluation posi- 
tions for the shear within them. In Section 3.4 we explain 
how the computed shear values are converted (after inte- 
gration along the lines of sight) to physical units, and the 
approximations we have adopted in the method. 

In Section 4.1 the Dyer-Roeder equation for the angu- 
lar diameter distances is introduced and Section 4.2 sum- 
marises the results derived in the Appendix for the beam 
equation generalised for all the cosmologies. This precedes 
the discussion (Section 4.3) of magnifications obtainable in 
inhomogeneous universes with different degrees of smooth- 
ness, and highlights the significance of differences between 
researchers who adopt different approaches to this subject. 

Section 5 formally presents the weak lensing results. 
We first (Section 5.1) attempt to see if the onset of struc- 
ture formation can be seen from the shear data, and briefly 
comment on the expected behaviour of the developing shear. 
Secondly (Section 5.2), the results for magnification, con- 
vergence, shear and ellipticity are presented for all the cos- 
mologies, compared and contrasted. We consider the impact 
of the smoothness parameter on the results. 

It may be possible to determine the density parameter, 
Oo, from the probability distribution and the skewness in 
the distribution for the convergence, both of which may be 
measurable observationally. In Section 6.1 we present our 
results for these quantities, and compare them with others. 
The application to the determination of the cosmological pa- 
rameters from Type la Supernovae is mentioned in Section 
6.2. 

In Section 7 the weak lensing statistics detailed in 
Section 5 are summarised, together with the results of the 
non-Gaussianity in the convergence. We also compare and 
contrast our results with other authors. Of particular signif- 
icance are the common use of two-dimensional approaches 
by others, and the use of either point mass particles or parti- 
cles with small softening scales, which introduce high values 
of magnification. These authors often use the empty cone 



approximation, rather than the full beam approximation, 
making comparisons difficult. This latter point is discussed. 

In the Appendix we state the generalised beam equa- 
tion, and derive the equations necessary for the numerical 
determination of the angular diameter distances for all the 
cosmologies. 

1.2 Previous work 

Numerous methods have been employed to study weak lens- 
ing, and throughout this paper we will make comparisons 
with previous work by other authors. We summarise here 
the methods which have been used by others. 

Jaroszyhski et al. (1990) use a 'ray-tracing' method to 
evaluate the matter column density in a matrix of 128 3 
pixels for each of their lens-planes. The boxes were gen- 
erated using a particle-mesh (PM) code in the standard 
Cold Dark Matter (SCDM) cosmology, and were of side- 
dimension 128/i _1 Mpc, where h is the Hubble constant in 
units of 100 km s _1 Mpc -1 . By making use of the assumed 
periodicity in the particle distribution orthogonal to the line 
of sight, they translate the planes for each ray, so that it be- 
comes centralised within the plane of one full period in ex- 
tent. This removes any bias acting on the ray when the shear 
is computed. Instead of calculating the effect of every parti- 
cle on the rays, the pixel column densities in the single pe- 
riod plane are used, and they assume that the matter in each 
of the pixels resides at the centre point of each pixel. They 
calculate the two two-dimensional components of the shear 
(see Section 2 for the definition of shear) as ratios of the 
mean convergence of the beam, which they obtain from the 
mean column density. However, they have not employed the 
net zero mean density requirement in the planes, (described 
in detail by Couchman, Barber and Thomas, 1999), which 
ensures that deflections and shear can only occur when there 
are departures from homogeneity. Also, the matter in the 
pixel through which the ray is located is excluded. To fol- 
low the shearing across subsequent planes they recursively 
generate the developing Jacobian matrix for each ray, in ac- 
cordance with the multiple lens-plane theory (see Section 
2). 

Wambsganss, Cen and Ostriker (1998) also use the 'ray- 
tracing' method in cosmological iV-body simulations. The 
method is applied to PM simulations of the SCDM cos- 
mology, in which a convolution method is used to combine 
large-scale boxes of 400/i _1 Mpc and resolution 0.8/i _1 Mpc, 
with small-scale boxes of 5/i _1 Mpc and the higher resolu- 
tion of 10/i _1 kpc. They randomly orient each simulation 
box, and project the matter contained within each onto a 
plane divided into pixels. They choose the central 8/i -1 Mpc 
x 8h~ Mpc region through which to shoot rays, but account 
for the deflections of the rays in terms of all the matter in the 
plane of 80/i _1 Mpc x 80/i _1 Mpc. However, to speed up the 
computation, a hierarchical tree code in two dimensions is 
used to collect together those lenses (pixels) far away, whilst 
treating nearby lenses individually. The code assumes that 
all the matter in a pixel is located at its centre of mass. The 
matter in each pixel, which measures 10ft _1 kpc x 10/i _1 kpc, 
is assumed to be uniformly spread. By using the multiple 
lens-plane theory, they show both the differential magnifi- 
cation probability distribution, and the integrated one for 
100 different source positions at redshift z s = 3.0. One ad- 
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vantage of this type of ray-tracing procedure is its ability 
to indicate the possibility of multiple imaging, where differ- 
ent rays in the image plane can be traced back to the same 
pixel in the source plane, and they are able to compute the 
statistics of angular separations for multiple images. 

Marri and Ferrara (1998) use lens-planes up to red- 
shifts of z = 10 for mass distributions determined by 
the Press- Schechter formalism, and treat the particles as 
point-like masses with no softening. They apply their ray- 
tracing method to three cosmologies with (Qm, ^a, Cl v ) = 
(1, 0, 0), (0.4, 0.6, 0), and (0.7, 0, 0.3). (Q M , Q\ and Sl v 
represent the density parameters for matter, vacuum en- 
ergy and the hot Dark Matter component, respectively.) 
The maximum number of lenses in a single plane is approx- 
imately 600, each having the appropriate computed mass 
value, and they follow 1.85 x 10 7 uniformly distributed rays 
within a solid angle of 2.8 x 10~ 6 sr, corresponding to a 
420" x 420" field. The final impact parameters of the rays are 
collected in an orthogonal grid of 300 2 pixels in the source 
plane. Because of the use of point masses, their method pro- 
duces very high magnification values, greater than 30 for the 
Einstein-de Sitter cosmology. They have also chosen to use 
a smoothness parameter a = in the redshift-angular diam- 
eter distance relation (described in Section 4) which depicts 
an entirely clumpy universe. 

Jain, Seljak and White (1998, 1999) have made use 
of iV-body simulations generated using a parallel adaptive 
particle-particle, particle-mesh (AP 3 M) code. The cosmolo- 
gies simulated are the same as the ones reported on here, but 
they have 256 3 particles and comoving box side dimensions 
of 85/t _1 Mpc for their Einstein-de Sitter cosmologies, and 
141/i~ Mpc for the open and cosmological constant mod- 
els. For each plane, the projected density, together with 
the appropriate redshift-dependent factors, is Fourier trans- 
formed, using periodic boundary conditions, to obtain the 
shear in Fourier space. The two-dimensional shear matrix 
is then computed onto a grid in real space by the inverse 
Fourier transform. The size of the grid is chosen to be com- 
patible with the force-softening scale (the resolution) in the 
iV-body simulations, so that for sources at redshifts around 
1, the angular-scale size of the grid is less than, or of the 
same order as, the angular-scale size of the force-softening 
at these redshifts. Perturbations on the photon trajectories 
are computed and the shear matrix interpolated to the pho- 
ton positions enabling the Jacobian matrices to be computed 
by recursion. In view of the use of a fine grid for the shear 
and deflection angle computations, they are able to analyse 
their data on different scales, and are thus able to determine 
the power spectrum in both the shear and the convergence. 
This approach is very different from the approach we fol- 
low which makes use of variable softening for the particles 
depending on their particular environment, and which en- 
sures that most 'rays' would pass entirely through regions 
of smoothed density, thereby requiring the application of 
the full beam approximation. Jain, Seljak and White (1998, 
1999) assume point particles interpolated onto their grid, 
and use the empty cone approximation. 

Hamana, Martel and Futamase (1999) study weak 
lensing in iV-body simulations of three cosmologies with 
(fi , A ) = (1, 0), (0.3, 0.7), and (0.3, 0). (We now use 
fio and Ao to represent the present matter density and vac- 
uum energy density parameters respectively.) These were 



generated by a P M algorithm, using Fourier techniques 
on a 128 3 lattice. The comoving simulation box sizes were 
128Mpc, and the particles were given comoving softenings 
of 300kpc. Three simulations were performed for each cos- 
mological model, and boxes from the different simulations 
combined to limit correlations in large-scale structure be- 
tween adjacent boxes. Each of the mass distributions in the 
simulation boxes were than projected onto planes at the box 
redshift, and Poisson's equation solved numerically on each 
of them. This was done by first evaluating the surface den- 
sity onto a 512 x 512 grid, based on the particle positions, 
and then inverting Poisson's equation using a Fast Fourier 
Transform (FFT) method. As the multiple lens-plane theory 
was to be used for more than 10 7 rays passing through each 
plane, the recursion algorithm was simplified by assuming 
small deflections for the rays. This also meant that the rays 
could be considered to pass through the grid points, and 
travel, effectively, in straight lines through the entire dis- 
tance from the observer to the source plane. We have also 
made this approximation in our own method. 

An alternative to the conventional form of 'ray-tracing' 
was introduced by Refsdal (1970), who used 'ray-tracing' 
with calculations of the differential deflections of light rays 
around a central ray to determine the distribution of mag- 
nifications. 

Fluke, Webster and Mortlock (1999) and Fluke, Web- 
ster and Mortlock (2000) have developed this idea further as 
the 'ray-bundle' method. The principle is to trace the pas- 
sage of a discrete bundle of light rays as it passes through 
the deflection planes. The advantage of the method is that it 
provides a direct comparison between the shape and size of 
the bundle at the observer and at the source plane, so that 
the magnification, ellipticity and rotation can be determined 
straightforwardly. The authors have selected a range of pop- 
ular cosmologies, including the ones reported on here, and 
have produced their own iV-body simulation data-sets using 
a P M algorithm. For each cosmology they have produced a 
number of independent simulations to enable them to ran- 
domly choose a simulation box for a particular epoch from 
any of the realisations. In this way, correlations of large-scale 
structure between adjacent boxes are avoided. They have 
run their simulations with 64 3 dark matter particles, with 
box sizes ranging from 80/i _1 Mpc to 164.3/i _1 Mpc. Having 
randomly selected the boxes for a given cosmology, the par- 
ticle mass distributions are then projected onto planes at the 
redshifts of the boxes. The authors construct approximately 
50,000 bundles, each comprising 8 rays, to shoot in random 
directions through the planes, from the observer's location 
at z — 0. The shooting area was limited to 50" x 50" to avoid 
edge effects of the planes, and only matter within a single 
period in the transverse direction is included in the planes. 
In addition, a radius, typically 15ft _1 Mpc centred around 
each bundle is chosen for the extent of the matter to be in- 
cluded in the deflection angle computations. The projected 
masses are also considered as point particles, so that very 
high magnification values can be achieved in principle; how- 
ever, the authors do not include bundles which pass within 
\/2 of the Einstein radius of any particle. Because of their 
use of point masses, they use the empty cone approxima- 
tion in the determination of the statistical distributions of 
magnifications for the different cosmological models. 

Premadi, Martel and Matzner (1998a) have introduced 
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individual galaxies into the computational volume, match- 
ing the 2-point correlation function for galaxies. They also 
assign morphological types to the galaxies according to the 
individual environment, and apply a particular surface den- 
sity profile for each. Five different sets of initial conditions 
were used for the simulations, so that the individual plane 
projections can be selected at random from any set. TV- 
body simulations were produced for three cosmologies with 
(fi 0) Ao) = (1, 0), (0.2, 0), and (0.2, 0.8) in boxes of comoving 
side-dimension 128/i _1 Mpc. They solve the two-dimensional 
Poisson equation on a grid and invert the equation using a 
FFT method to obtain the first and second derivatives of 
the gravitational potential on each plane. They also cor- 
rectly ensure that the mean surface density in each lens- 
plane vanishes, so that a good interpretation of the effects 
of the background matter is made. Their method uses beams 
of light, each comprising 65 rays arranged in two concentric 
rings of 32 rays each, plus a central ray. The multiple lens- 
plane theory then enables the distributions of cumulative 
magnifications to be obtained. 

Tomita (1998a and b) also uses a ray bundle method, 
but by evaluating the gravitational potential at approxi- 
mately 3000 locations between the observer and sources at 
redshift 5, he is able to compute the weak lensing statis- 
tics without using the multiple lens-plane theory. He has 
used TV-body simulations produced using a tree code with 
32 3 particles in four different cosmologies with (Qo, Ao) = 
(1, 0), (0.2, 0.8), (0.2, 0), and (0.4, 0). In all of these, except 
the SCDM cosmology (fio = 1, Ao = 0), the particles are 
treated as galaxy-size objects. In the SCDM cosmology, 20% 
of the particles are treated thus, whilst 80% are given soften- 
ing radii of 100/i~ kpc (model a), 40ft _1 kpc (model b), and 
20h~ kpc (model c) . Ray bundles are formed from 5x5 rays 
arranged in a square, with separation angles (fixed for each 
run) at between 2 arcsec and 1 degree. At each of the 3000 
positions between the observer and source, the potential is 
evaluated at each ray position by translating each simula- 
tion cube (using the periodic properties) so that the centre of 
each bundle is located in the centre of each cube. In this way, 
all the matter within a full period in extent contributes to 
the calculation of the potential (although no account is taken 
of matter beyond one period transverse to the line of sight). 
To avoid spurious values of the potential arising from masses 
close to evaluation positions, an average of the potential is 
taken by integrating it analytically over the interval between 
adjacent evaluation positions. The light propagation is then 
determined by solving the null-geodesic equations, and the 
required statistics constructed from shooting 1000 bundles 
through the flat cosmologies, and 200 bundles through the 
open cosmologies. 



2 THE PROPAGATION OF LIGHT 

In the case of multiple deflections of light by a series of pro- 
jected lens-planes, the Jacobian matrix develops in accor- 
dance with the multiple lens-plane theory, which has been 
summarised by Schneider, Ehlers and Falco (1992). Its form 
at redshift z = enables all the lensed properties of light 
from distant sources to be determined. 

We follow, in outline, their description. For a system of 
TV lenses, the basic lens equation for a single lens, relating 



the angular position of the source to that of the image, may 
be generalised straightforwardly. If the position vector of 
the source in the plane perpendicular to the line of sight at 
the source position (the source plane) is rj, and the position 
vectors of the image positions in the various, TV, deflection 
planes are £i, where i — 1, N , then the lens equation may 
be written as 

JV 

where Di is the angular diameter distance to the ith lens, 
Di s is that from the ith lens to the source, and &i is the 
deflection angle at the ith lens. To make equation |l| dimen- 
sionless, put 



Xi = ti/Di, 

and for the source, 

y = xjv+i = 7]/D s = tj/Dn+i. 



(2) 



(3) 



We also use the individual reduced deflection angles defined 
by 



OLi — -J—OLi. 

Us 

Then the displacement in the jih lens plane is 



j-i 



i-i 

, , \ ~* D s D 

-a i (fl«)=x 1 -^-- 



Xj = xi— ^Y~ai{DiKi) = xi- 



•j 



Xl 



3-1 



f3ijOn(DiXi), 



where 

Pij - D~ s Dj ■ 



(4) 



«(DiXi),(5) 



(6) 



(7) 



Then the full form of the ray-tracing equation is equation ^ 
evaluated at the source plane: 



Xl 



N 

E 

1=1 



a,. 



(8) 



(The /3ij factor disappears, because /3j s = 1 from its defi- 
nition.) The mapping of the source onto the image is given 
by the Jacobian matrix, which relates small changes in the 
source to corresponding small changes in the image seen on 
the first lens-plane: 

dy 



A 



dxi 



(9) 



Similarly, we may define Jacobian matrices appropriately at 
each of the lens-planes: 

9xi 



A, 



<9xi 



(10) 



We now define the derivative of the reduced deflection angle 
for the ith lens by 

don 



11 



<9xi 



(11) 



By defining the 'effective lensing potential' in terms of the 
angular position, Oi for the ith lens as an integral along the 
line of sight direction, X3, 
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D d 



D d D„ 



(f>(D d 0i,x 3 )dx3, 



(12) 



in which c is the velocity ol light, Ui can be shown to be 
equivalent to the matrix of the second derivatives of the 
effective lensing potential for the ith lens, and is therefore 
related to the second derivative of the gravitational potential 
(the shear): 



Ui = 



4>ii V>i2 

V4l V>22 



(13) 



where the superscripts, i, denote the deflection plane index, 
and where we have written 



d 2 ^{6) 
d 2 i>{0) 



d 2 ^(6) 



and tp22 



hd9 2 ' 



(14) 



(15) 



de 2 de 1 

in which the suffixes in the denominators refer to the coor- 
dinate directions. Then the ray-tracing equation (fel) gives 



A 



dy 

C*Xi 



dog 
^ <9xi 

\ -» da. i 9xi 



A' 



j=l i=l 



(16) 



where T is the identity matrix. 

Thus the final Jacobian matrix can be evaluated at z = 
0, since the individual matrices can be obtained by recursion. 
Using equation H, they are just 



for the jth lens, and 
Ai =T. 



(17) 



(18) 



In our approach, the second derivatives of the two- 
dimensional effective lensing potentials required for each 
deflection location, are obtained by integration of the com- 
puted three-dimensional shear values, i.e., the second deriva- 
tives of the peculiar gravitational potential. It is necessary 
to work with the peculiar gravitational potential, because 
the shearing of light arises from deviations from homogene- 
ity; in a pure Robertson- Walker universe we would expect 
no deviations. Couchman, Barber and Thomas (1999) de- 
rive the expression for the peculiar gravitational potential, 
(f>, in terms of the gravitational potential, $, and the mean 
density, p: 



= ^ -^j^Gapx 2 ; 



(19) 

G is the universal gravitational constant, x is the position 
vector, and a is the expansion factor for the universe (so 
that ax is the comoving position vector). This result corre- 
sponds to a system with zero net mass on large scales and 
immediately gives 



dxi 
and 
d 2 



dx i 



- 4 / 3 nGa pxi 



d 2 $ 



dxidxj dxidxj 



4/ 3 nGa pSij. 



(20) 



(21) 



Equation |2l] shows how the real result for the shear, 
d 2 <p/dxidxj, based on the peculiar gravitational potential, is 
related to the value of d 2 § / dxidxj through the subtraction 
of the term in the mean density. 

Equation |2l] can now be evaluated explicitly for the 
three-dimensional shear. The integral solution of Poisson's 
equation is well-known, and the solution can be easily dif- 
ferentiated twice to give 

<9 2 <XR) _ 



dxidxj 
f I 

G 



p(R') 



R — R' 



3p(B!)(x i -x' i )(x j -x' j ) 



| R-R | 5 

4 / 3 TvGa 2 p5i 



d 3 R' 



(22) 

(We have introduced R and R' for the evaluation posi- 
tion for the shear, and the matter positions respectively. In 
practice, of course, the triple integral over all space would 
be replaced by a summation.) The two-dimensional second 
derivatives of the effective lensing potentials required for the 
Jacobian matrices then follow immediately from equation |l2| 
(using spatial rather than angular coordinates): 



D d D d 



d 2 </>(x 3 ) 
dxidxj 



dx 3 



(23) 



Dd, D ds , and D s are the angular diameter distances from 
the observer to the lens, the lens to the source, and the 
observer to the source, respectively. The integration is along 
the coordinate direction, x 3 , and the subscripts i and j now 
refer to any of the three coordinate directions. This equation 
applies quite generally for any deflection location, so the 
deflection plane index has been dropped for clarity. 

We now summarise the main equations we shall need for 
weak lensing, which are obtainable from the final Jacobian 
matrix. The final emergent magnification, fi, may be com- 
puted after passage through an entire box or set of boxes, 
and is 



p = (det A) 1 ■ 

The convergence, k, is 

K = gW'll + ^22), 



(24) 



(25) 



and is therefore obtainable from the diagonal elements of 
the Jacobian matrix. 

The two-dimensional shear, 7, in each line of sight, is 
given by 



7 = 7(^11 ^ ^22) + ipi 



(26) 



(We may take 1^12 = ^21, because we are dealing with a 
weak shear field which is smoothed by the variable particle 
softening, ensuring that the gravitational potential and its 
derivatives are well-behaved continuous functions.) 
From equation and these definitions, 

p = (1 - Vn - ^22 + ^11^22 - ^i2)~\ (27) 



(28) 



^ (1-k) 2 -7 2 ' 

In the presence of convergence and shear, a circular source 
becomes elliptical in shape, and the ellipticity, e, defined in 
terms of the ratio of the minor and major axes, becomes 
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(29) 

1 - k + 7 
which reduces to 

e~2 7 (l + K -7)+0( K 3 , 7 3 ) (30) 

in weak lensing. 

In our method, the evaluation of the second deriva- 
tives of the two-dimensional effective lensing potentials is ob- 
tained from integration of the computed three-dimensional 
shear values at a large number of evaluation positions along 
lines of sight. The multiple lens-plane theory then enables 
distributions of the magnification, ellipticity, convergence 
and shear at redshift z = to be computed for light rays 
traversing the set of linked simulation boxes starting from 
the chosen source redshift. The ability to apply the appro- 
priate angular diameter distances at every evaluation posi- 
tion avoids the introduction of errors associated with planar 
methods, and also allows the possibility of choosing source 
positions within a simulation box if necessary. 

3 PROCEDURE 

3.1 The three-dimensional shear algorithm 

Couchman et al. (1998) describe in detail the algorithm for 
the computation of the elements of the matrix of second 
derivatives of the gravitational potential. The algorithm is 
based on the standard P 3 M method (see Hockney and East- 
wood, 1988), and uses a Fast Fourier Transform convolution 
method. It computes all of the six independent shear compo- 
nent values at each of a large number of selected evaluation 
positions within a three-dimensional TV-body particle sim- 
ulation box. It has a computational cost of order Nlog 2 N, 
where N is the number of particles in the simulation volume, 
and for ensembles of particles, used in typical iV-body sim- 
ulations, the rms errors in the computed shear component 
values are typically ~ 0.3%. In addition, the shear algorithm 
has the following features. 

A) The algorithm uses variable softening designed to 
distribute the mass of each particle within a radial profile 
depending on its specific environment. By virtue of this fa- 
cility, we have been able to choose the softening such that 
light rays feel the existence of a smooth mass distribution. 
Each particle may be assigned its own softening-scale pa- 
rameter, depending on the particle number-density in its 
environment. In this way, it can be used to minimise the ef- 
fects of isolated single particles, whilst the smoothed denser 
regions are able to represent the form of the large-scale struc- 
ture. The parameter we have chosen to delineate the soft- 
ening scale for each particle is proportional to I, where 21 is 
the radial distance to the particle's 32nd nearest neighbour. 
The value of I has been evaluated for every particle by using 
a smoothed particle hydrodynamics (SPH) density program 
for each simulation box, and is read in by the shear code 
along with the particle position coordinates. 

The maximum softening is allowed to be of the order of 
the mesh dimension for isolated particles, which is defined 
by the regular grid laid down to decompose the short- and 
long-range force calculations. In this way individual isolated 
particles are unable to strongly influence the computed shear 
values, in accordance with our need to study the broad prop- 
erties of the large-scale structure, rather than the effects of 



individual particles, which are not representative of physical 
objects. 

To determine realistic values for the minimum softening 
scale, we wanted to keep the incidence of strong lensing to 
a minimum, whilst at the same time allowing a physical in- 
terpretation. Using the dimensionless angular diameter dis- 
tances in terms of the present value of the Hubble parameter, 
H , i.e., r d (= Sa.D d ), r ds (= ^D ds ), and r s (= ^D s ), 
for the observer to the lens, the lens to the source, and the 
observer to the source, respectively, the Einstein radius, Re, 
becomes 

Re = 8.6 x 10- s ft _1 JV» (j^) 5 Mpc (31) 

for a cluster of N particles, each of mass 1.29 x 10 11 solar 
masses (see Section 3.2), where ft is the Hubble parameter 
in units of 100 km s _1 Mpc -1 . Substituting values for the 
angular diameter distance factors then gives a maximum 
value of Re = 2.84 x 10 -3 ft -1 iV2 Mpc for a source at red- 
shift z s = 1 in the SCDM cosmology. This occurs for a lens 
at redshift z d = 0.29. Thus, for a cluster of 1000 particles, 
Re = 0.089ftT 1 Mpc. For a source at redshift 3.6, the max- 
imum value of Re is 0.108ft. -1 Mpc and occurs for a lens at 
redshift 0.53. Thus by setting a working minimum value for 
the variable softening of 0.1ft -1 Mpc, we would rarely ex- 
pect to see strong lensing due to caustics in the simulations. 
In box units, the softening is 10 x (1 + z), where z is the 
box redshift. 

The corresponding values for Re in the other cosmolo- 
gies, for N = 1000, are as follows. 

For flo = 0.3, Ao = 0, the maximum value of Re is 
0.093ft - Mpc for z s — 1, and occurs for z d = 0.32. For 
z s = 3.6, the maximum value of Re is 0.115ft -1 Mpc, and 
occurs for z d — 0.58. 

For O,o = 0.3, Ao = 0.7, the maximum value of Re is 
0.104ft~ 1 Mpc for z s = 1, and occurs for z d = 0.36. For 
z 3 = 3.6, the maximum value of Re is 0.142ft -1 Mpc, and 
occurs for z d — 0.67. 

Consequently, we are able to ensure that the incidence 
of strong lensing is kept to very low levels in all cosmologies, 
and that the minimum softening is always greater than, or 
similar to, the maximum value of the Einstein radius for a 
cluster of 1000 particles. At the same time, the softening 
scale is approximately of galactic dimensions, giving a re- 
alistic interpretation to the choice. We have also accounted 
for the effects of the different numbers of particles per box 
in the different cosmological simulations. The variable soft- 
ening scale for each particle has been retained at the same 
level in all the cosmologies, so that the same mass value 
is contained within it. Thus, around small-particle clusters 
there would be differences in the softening from cosmology 
to cosmology, but around large clusters, where the shearing 
is likely to be most important, the differences would be tiny. 

Couchman et al. (1998) have investigated the sensitiv- 
ity of weak lensing results to the input softening, finding 
differences only in limited numbers of lines of sight at the 
high magnification end of the distributions. 

B) The shear algorithm works within three- 
dimensional simulation volumes, rather than on planar pro- 
jections of the particle distributions, so that angular diame- 
ter distances to every evaluation position can be applied. 
It has been shown (Couchman et al., 1998) that in spe- 
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cific circumstances, the results of two-dimensional planar 
approaches are equivalent to three-dimensional values inte- 
grated throughout the depth of a simulation box, provided 
the angular diameter distance is assumed constant through- 
out the depth. However, by ignoring the variation in the 
angular diameter distances throughout the box, errors up 
to a maximum of 9% can be reached at a redshift of z = 0.5 
for SCDM simulation cubes of comoving side 100h~ Mpc. 
(Errors can be larger than this at high and low redshift, 
but the angular diameter distance multiplying factor for the 
shear values is greatest here for sources we have chosen at a 
redshift of 4.) 

C) The shear algorithm automatically includes the con- 
tributions of the periodic images of the fundamental vol- 
ume, essentially creating a realisation extending to infinity. 
Couchman et al. (1998) showed that it is necessary to include 
the effects of matter well beyond the fundamental volume 
in general (but depending on the particular particle distri- 
bution), to achieve accurate values for the shear. Methods 
which make use of only the matter within the fundamen- 
tal volume may suffer from inadequate convergence to the 
limiting values. 

D) The method uses the peculiar gravitational poten- 
tial, <f>, through the subtraction of a term depending upon 
the mean density. Such an approach is equivalent to requir- 
ing that the net total mass in the system be set to zero, and 
ensures that we deal only with light ray deflections arising 
from departures from homogeneity. 

3.2 The Hydra iV-body simulations 

The three-dimensional shear code can be applied to any 
three-dimensional distribution of point masses confined 
within a cubic volume, and produces shear values as if the 
fundamental volume were repeated indefinitely to represent 
a three-dimensional periodic distribution of masses. Each 
particle may be assigned an individual mass, although in 
the tests and our application of the code, all the particles 
were assumed to be dark matter particles with the same 
mass. 

The code has been applied to the data bank of cos- 
mological iV-body simulations provided by the Hydra Con- 
sortium ( jittp: / /hydra. mcmaster.ca/hydra/index. html ) and 
produced using the 'Hydra' iV-body hydrodynamics code 
(Couchman, Thomas and Pearce, 1995). Simulations from 
four different cosmologies were used, which will be referred 
to as the SCDM, TCDM, OCDM and LCDM cosmologies. 
Each of the simulations used a Cold Dark Matter-like spec- 
trum, and the parameters used in the generation and specifi- 
cation of these cosmological simulations are listed in Table |l| 

£7o and Ao are the present-day values of the density pa- 
rameter and the vacuum energy density parameter respec- 
tively, so that the SCDM and TCDM cosmologies are repre- 
sentative of Einstein-de Sitter universes, whereas the OCDM 
cosmology represents a low density, open universe, and the 
LCDM a low density, but spatially-flat universe with a cos- 
mological constant. The power spectrum shape parameter, 
r, is set to 0.5 in the SCDM cosmology, but the empirical 
determination (Peacock and Dodds, 1994) of 0.25 for cluster 
scales has been used in the other cosmologies. In each case, 
the normalisation, as, on scales of 8/i _1 Mpc has been set to 



Cosmology 


n 


Ao 


r 


CT8 


No. of 


Box side 












pjirtriclGS 


( h — 1 Mnr^ 


SCDM 


1.0 


0.0 


0.50 


0.64 


128 3 


100 


TCDM 


1.0 


0.0 


0.25 


0.64 


128 3 


100 


OCDM 


0.3 


0.0 


0.25 


1.06 


86 3 


100 


LCDM 


0.3 


0.7 


0.25 


1.22 


86 3 


100 



Table 1. Parameters used in the generation of the four different 
cosmological simulations. 



reproduce the number density of clusters (Viana and Liddle, 
1996). 

In the SCDM and TCDM cosmologies, the number of 
particles is 128 3 , leading to individual dark matter par- 
ticle masses of M part = 1.29 x lO 11 /!" 1 solar masses. In 
the low density universes, the number of particles is 0.3 
times the number in the critical density universes, leading 
to the same individual particle masses. The simulation out- 
put times were chosen so that consecutive simulation boxes 
may be snugly abutted; the side-dimensions are 100/i -1 Mpc 
in every case. Consequently, there are different numbers of 
time-outputs to a given redshift value for the different cos- 
mologies. For a nominal source redshift of z s — 4 (which is 
the furthest extent of our weak lensing analysis), 33 simula- 
tion boxes were abutted to a redshift of 3.90 in the SCDM 
cosmology, 33 to a redshift of 3.93 in the TCDM cosmology, 
41 to a redshift of 4.00 in the OCDM cosmology, and 48 to 
a redshift of 3.57 in the LCDM cosmology. (The number of 
boxes has no bearing on the weak lensing statistics.) 

3.3 Orientation and lines of sight 

Each time-output in a given simulation run is generated us- 
ing the same initial conditions, so that a particular structure 
(although evolving) occurs at the same location in all the 
boxes, and is therefore repeated with the periodicity of the 
box. To avoid such obvious and unrealistic correlations, we 
have arbitrarily translated, rotated (by multiples of 90°) and 
reflected each time-slice, about each coordinate axis, before 
linking them together to form the continuous depiction of 
the universe back to the source redshift. 

To follow the behaviour of light rays from a distant 
source through the simulation boxes, and obtain distribu- 
tions of the properties at z = 0, we construct a regular rect- 
angular grid of directions through each box. Since there are 
likely to be only small deflections, and the point of interest 
is the statistics of output values, each light ray is considered 
to follow one of the lines defined by these directions through 
the boxes. The evaluation positions are specified along each 
of these lines of sight. 

We have selected 1000 evaluation positions on each of 
100 x 100 lines of sight, which is well matched to the min- 
imum variable softening, giving adequate sampling in the 
line of sight direction; the method has been tested using up 
to a total of 4 x 10 6 lines of sight, and we have found that 
whilst this smooths the distribution plots for the magnifica- 
tion and shear at the high magnification end and gives rise 
to higher maximum values of the magnification, the statisti- 
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cal widths of the plots are virtually unchanged. Since we are 
dealing with weak lensing effects and are interested only in 
the statistical distribution of values, these lines of sight ade- 
quately represent the trajectories of light rays through each 
simulation box. It is sufficient also to connect each 'ray' with 
the corresponding line of sight through subsequent boxes in 
order to obtain the required statistics of weak lensing. This 
is justified because of the random re-orientation of each box 
performed before the shear algorithm is applied. 

3.4 Conversion factors and approximations in the 
method 

The second derivatives of the two-dimensional effective lens- 
ing potentials are obtained from the three-dimensional sec- 
ond derivatives of the peculiar gravitational potential by in- 
tegration, in accordance with equation The integration of 
the three-dimensional shear values has been made in small 
steps (0.02 of the box depth) along each line of sight, en- 
abling the weak lensing properties to be determined from 
the Jacobian matrices and recorded at 50 evenly-spaced lo- 
cations along each line of sight in every simulation box. To 
evaluate the absolute second derivatives of the effective lens- 
ing potentials, the appropriate scaling factor is introduced, 
which applies to the simulation box dimensions. From equa- 
tion^, the factor B(l + z) 2 rdrds/r a can be extracted, where 
B = (c/J c fo)(2/c 2 )GM part x (comoving box depth)" 2 . For 
the simulation boxes used, which have comoving dimensions 
of 100/i _1 Mpc, B = 3.733 x 10" 9 , and the (l + z) 2 factor oc- 
curs to convert the comoving code units into physical units. 
(The step length for the integrations in each box is left as a 
code parameter, so that it may be varied at will, although 
in the analysis reported here, the step length was invariably 
0.02 x the box depth.) 

The procedure to obtain the two-dimensional second 
derivatives of the effective lensing potentials involves a num- 
ber of approximations. 

First, it is assumed that the angular diameter distances 
required for equation ^ vary linearly throughout each step 
length (0.02 x the box depth). They are, however, evaluated 
exactly at the 50 step positions through each box. 

Second, in our approach we assume that, although each 
simulation box is generated as a single simulation output- 
time representation, the angular diameter distances vary 
throughout the depth of each box, as they would in the real 
universe. 

Third, as mentioned in the previous section, we assume 
only weak deflections for light rays, so that they are assumed 
to follow the straight lines of sight represented by the grid 
points of the evaluation positions. Since we are interested 
only in the statistics of lensing within large-scale structure 
simulations, this is a perfectly acceptable practice. 

Finally, we have made use of a weak lensing approxima- 
tion for the computation of the intermediate Jacobian matri- 
ces. The full form for the final Jacobian matrix (equation |l6|) 
is computed by recursion, by including the intermediate ma- 
trices given by equations ^ and Expanding equation ^ 
fully involves second and higher order terms, arising from 
the cross multiplication of the intermediate Jacobian matri- 
ces, the distance factors represented by and the matrices 
Ui. For weak lensing in which all the tpij are much less than 
unity, the full form for .Atotal simplifies considerably to 



-4totai — 1 — y ' Ui . 



(32) 



This is the form used; strong lensing events will still be 
recorded as such, but using this approximation, the com- 
ponent values in .Atotal will n °t be accurate in strong lens- 
ing cases. However, the incidence of strong lensing events is 
likely to be very small, because of our choice of the softening 
scale (see Section 4.3). This being the case, any strong lens- 
ing events will not adversely affect the weak lensing statistics 
determined in our analysis. 



4 ANGULAR DIAMETER DISTANCES IN 
INHOMOGENEOUS UNIVERSES 

4.1 The Dyer-Roeder equation 

Our three-dimensional approach allows the use of the appro- 
priate angular diameter distances at every single evaluation 
position. This is not possible in two-dimensional approaches, 
where it is assumed that all the lensing mass in a box is pro- 
jected onto a plane at a single angular diameter distance. 

Since the angular diameter distances depend very much 
on the distribution of matter and the particular cosmology, 
it is therefore necessary to have available appropriate values 
for the angular diameter distances for the particular distri- 
bution of matter in the simulation data-set being investi- 
gated. 

By considering the universe to be populated by ran- 
domly distributed matter inhomogeneities, but resembling 
the Robertson- Walker, Friedmann-Lemaitre model on large 
scales, a second order differential equation is obtained for 
the angular diameter distance, D, in terms of the density 
parameter, fio, for the universe, the vacuum energy density 
parameter, Ao, and the redshift, z, of the source. Dyer and 
Roeder (1973) made assumptions about the type of matter 
distribution to obtain a more practical equation for Ao = 
cosmologies. They assumed that a mass fraction, d, (called 
the smoothness parameter), of matter in the universe is 
smoothly distributed, and that the fraction (1 — a) is bound 
into clumps. Then the equation for the angular diameter 
distance (with Ao = 0) is: 



(* + l)[fio* + l] 



(£D 



+ 



dD 

dz 



+ 



-afi + 



1 + zY> 



L> = 0, 



(33) 



in which a is the optical scalar for the shear, introduced by 
the matter distribution around the beam. 

In order to apply equation p3|, Dyer and Roeder (1973) 
considered the following scenarios. They considered a uni- 
verse in which all the matter is bound into clumps, so that 
a = 0, and in which the light beam passes far away from the 
clumps. This is described as light propagating through an 
'empty cone,' and gives rise to maximal divergence of the 
beam. The opposite scenario has 5=1, i.e., an entirely 
smooth universe. Here the smooth matter distribution is 
present within the beam, giving a 'full cone,' or 'filled beam' 
approximation. 

We have considered whether the shear along individual 
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Figure 1. The angular diameter distance, for the different 
cosmologies, assuming a source redshift of z s = 3.6, and a = 1. 



Figure 2. The angular diameter distance multiplying factor, 
r d r ds/ r s in the different cosmologies, assuming a source redshift 
of z s = 3.6, and a = 1. 



lines of sight is able to significantly affect the chosen val- 
ues for the angular diameter distances. Our work has been 
conducted using cosmological simulations in which the dis- 
tributions of matter are smooth, and we show in Section 5.2 
that the effects are negligible. With a ~ 0, therefore, equa- 
tion [33] immediately reduces to the well-known Dyer-Roeder 
equation, 



+ 1) [n z + 1] ^ + 



2 



+ 3 



dD 

dz 



which can be solved analytically for fio 
arbitrary a. 



1, A 



(34) 
0, and 



4.2 Generalisation of the Dyer-Roeder equation 

Starting from the generalised beam equation, quoted by Lin- 
der (1998a and b), we have generalised the form of the Dyer- 
Roeder equation to apply to all the cosmologies we have sim- 
ulated. This was necessary because solutions for the angular 
diameter distances were required in the LCDM cosmology 
containing a vacuum energy density. The procedure to gen- 
eralise the Dyer-Roeder equation is described fully in the 
appendix, and solutions of the final equation were obtained 
numerically. Figure [j] shows the result of solving the gener- 
alised equation, with a = 1, in the different cosmologies, for 
a source redshift of z s = 3.6, and Figure ^ shows the values 
of the angular diameter distance multiplying factor, rdrds/rs 
(which we now denote by R), also for a = 1. It is clear from 
this plot that the angular diameter distance multiplying fac- 
tor is considerably higher in the LCDM cosmology than the 
other cosmologies, and we shall comment further on this in 
regard to the weak lensing statistics in the discussion of our 
results in Section 7. When a — 0, the values are lower than 
for a = 1; we have tabulated the ratios R(a = 1)/R(a — 0) 
for the different cosmologies in Table |^. 



Cosmology 


Minimum a 


R(a = 1)/R[a(z = 0)] 


SCDM 


0.83 


1.0408 


TCDM 


0.88 


1.0204 


OCDM 


0.80 


1.0241 


LCDM 


0.82 


1.0414 



Table 2. The values of a at 2 = and the ratios of the R factors 
at the peaks of the curves in the different cosmologies, assuming 
a source redshift of 3.6. 



4.3 Magnification in inhomogeneous universes 

Later, in Section 7, we make comparisons of our results with 
other workers, who may use either the full beam or empty 
beam approaches for the propagation of light. We find it dif- 
ficult to make meaningful comparisons with results obtained 
using the empty beam approach, because the magnification 
distributions, for example, may be quite different depending 
on the approach used. We therefore now consider the effects 
of inhomogeneities, in the different approaches, which are 
described loosely in terms of Dyer and Roeder's smooth- 
ness parameter, a. We follow the line of reasoning given by 
Schneider, Ehlers and Falco (1992). 

Consider our inhomogeneous universe to be on-average, 
i.e., on large scales, homogeneous and isotropic, so that the 
average flux from a source at redshift z and luminosity L 
will equal the flux, Sfl, observed in a smooth Friedmann- 
Lemaitre universe without local inhomogeneities: 



4n[D L (z)Y 



(35) 



Dl{z) is the luminosity distance in the smooth Friedmann- 
Lemaitre model, and is the mean of the luminosity dis- 
tance values in an inhomogeneous universe including the 
constraints of flux conservation. Thus Dl{z) can be related 
to the Dyer-Roeder angular diameter distance, D s , in an 
entirely smooth universe (with 5 = 1): 
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D L {z) = {l + z) 2 D s {z;a = l). 



(36) 



Now the magnification, fi, is just the ratio of the flux actu- 
ally observed in the image of a source and the flux which 
the same source would produce if observed through an 
empty cone without deflection. Then equation straight- 
away gives, for the mean magnification in terms of the ap- 
propriate Dyer-Roeder angular diameter distances, 



0*) 



D s (z; a) 



D 3 (z;d = 1) 



(37) 



Clearly, the magnification values derived in this way depend 
on the approximation used, and specifically the value of a. 
For example, rays passing close to clumps or through high- 
density regions will result in magnification in any approx- 
imation. If the empty cone approximation is used, then /j, 
will be greater than 1, and if the full cone approximation is 
used, then /i will be greater than the mean magnification. 
From equation |37| it follows immediately that the mean mag- 
nifications, (fi e ) and (/x/), in the empty cone and full beam 
approximations respectively are 



(fie) 



and 



D 3 {z;a = G) 



D 3 (z;d = l) 



(38) 



(W) = !• (39) 

(Me) > 1 because D s (z;a — 0) > D s (z;d — 1) at all red- 
shifts in all of our cosmologies. This is an important result 
for advocates of the empty beam approximation, particu- 
larly those working with distributions of point mass parti- 
cles, because by evaluating numerically the angular diam- 
eter distance factors in the different cosmologies, (see the 
Appendix), it immediately follows that 



(At e )(OCDM) < ( Mc )(LCDM) < (/x e )(SCDM). 



(40) 



Pei (1993) succeeded in calculating the statistical proper- 
ties of the magnifications due to a random distribution of 
point mass lenses using the assumption that the total mag- 
nification is a result of multiplication of the magnifications 
produced at each redshift interval. He found that the mean 
magnification as a function of redshift was exponential in 
terms of the optical depth, t(z): 



{He(z)) = exp[2r(z)]. 



(41) 



(The optical depth is the fraction of the sky covered by cir- 
cles of Einstein radii between the observer and the specified 
redshift, and is therefore dependent on the cosmology or 
the distribution and density of matter. Pei (1993) gives an 
expression for the optical depth at redshift z in terms of 
the Dyer-Roeder smoothness parameter and the cosmologi- 
cal density parameter for the point mass lenses in the SCDM 
cosmology.) Although these are useful results, interpretation 
of the magnification distributions in the empty cone approx- 
imation for different cosmologies is often complicated by the 
high magnification tails in the distributions which arise from 
light rays passing close to point mass particles. Rauch (1991) 
set up a random distribution of point masses and performed 
Monte Carlo simulations to calculate the resulting amplifica- 
tion probability distributions. He then fitted the distribution 
by an analytical expression for the probability which was 
given only in terms of the mean magnification. Thus, for a 



£ 
3. 




Figure 3. The minimum magnification values versus redshift for 
the filled beam approximation in three different cosmologies. 



given mean magnification, the distribution curves would be 
identical, and almost certainly somewhat unrealistic. How- 
ever, since the mean magnification for sources at a given 
redshift is dependent on the cosmology, the actual distri- 
bution curves are able, in principle, to distinguish between 
cosmologies. 

The minimum magnification in the empty cone approx- 
imation is 



Me,mm 



1. 



(42) 



Thus rays passing through voids will have = /i e ,mm = 1 
in the empty cone approximation, since the rays will be far 
from all concentrations of matter, and will satisfy the empty 
cone conditions. 

In the full cone approximation however, the magnifi- 
cation for rays passing through voids will be less than or 
equal to 1 because the rays will suffer divergence, and the 
minimum value will be 



D{z;a = 1) 
D(z;d = 0) 



< 1. 



(43) 



From this equation, we can again make comparisons 
amongst our cosmologies. Figure plots the value of /j/ !m i n 
from this equation versus redshift for the different cosmolo- 
gies, and immediately we can see that 

Li f ,min(SCDM) < w , min (LCDM) < w , min (OCT>M) (44) 

for all redshifts. 

There are two points to note about this result. First, 
the minimum magnification values will only rarely (if ever) 
be seen, because this would require the rays to pass entirely 
through the most underdense regions. For this reason the 
values may only be treated as lower bounds for the com- 
puted values. Comparisons between the computed minima 
and these analytic values are shown in Table || As required, 
the computed values are consistently greater than the the- 
oretical minima. Second, the values say nothing about the 
distribution of magnifications in any cosmology; it is not nec- 
essary, for example, for the cosmology producing the lowest 
minimum magnification to have the broadest range for the 
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Redshift 


//f (analytic) 


Uf fcomDutcd) 


SCDM/TCDM 






0.5 


0.9245 


0.9547 


1.0 


0.7892 


0.9030 


2.0 


0.5656 


0.8325 


4.0 


0.3201 


0.7486 


OCDM 






0.5 


0.9734 


0.9801 


1.0 


0.9147 


0.9510 


2.0 


0.7862 


0.8888 


4.0 


0.4192 


0.7922 


LCDM 






0.5 


0.9662 


0.9766 


1.0 


0.9264 


0.8791 


2.0 


0.6854 


0.8253 


4.0 


0.5802 


0.7145 



Table 3. The analytical minimum magnification values in the 
filled beam approximation for the different cosmologies, and the 
computed values from the simulations. It is necessary only that 
the computed values should be larger than the minima, since the 
theoretical minima will not, in practice, be seen. 



probability distribution. Consequently, the results reported 
in Section 5 require only that the minimum magnifications 
satisfy the minimum theoretical bounds, and may not re- 
late in any specific way to results of other workers using 
point mass particles in an empty cone scenario. This point 
is discussed further in Section 7. 



5 WEAK LENSING RESULTS IN THE 
DIFFERENT COSMOLOGIES 

5.1 The formation of structure 

The formation of structure occurs at different rates in the 
different cosmologies. Richstone, Loeb and Turner (1992), 
for example, considered the spherical collapse of density per- 
turbations, starting from an initial Gaussian distribution, 
and found that the rate of cluster formation as a function of 
redshift depended crucially on the value of Slo ■ This has been 
confirmed by, for example, Bartelmann, Ehlers and Schnei- 
der (1993). Lacey and Cole (1993), starting from the basic 
Press-Schechter formula, derived an equation for the merger 
rates of virialised halos in hierarchical models, which again 
showed the rates to be crucially dependent on fio- Later, 
Lacey and Cole (1994) compared their analytical results 
with the merger rates seen in iV-body particle simulations 
for the SCDM cosmology, and found good agreement for this 
cosmology. In addition, their analytical result was applicable 
to arbitrary values of f2o, and more general power spectra. 

Peebles (1993) summarises the evolution of structure in 
the Press-Schechter approximation, which provides the num- 
ber density for collapsed objects by mass scale. This can be 
evaluated in terms of the rms mass fluctuation, ay, at a fixed 
comoving scale, and so is redshift dependent. For Qo close 
to 1, ov oc (1 + z) -1 , and determines the evolution of the 
comoving number density of clusters. On this model, half 
of the present-day number of clusters would have formed 



later than z ~ 0.1, and 90% would have formed later than 
z ~ 0.3. This rapid evolution at late times in an Einstein- 
de Sitter universe is seen also in iV-body simulations. In 
low density universes the time evolution of oy is slower at 
low redshift, and this reduces the predicted rate of cluster 
formation at late times. The Press- Shechter approximation 
also underestimates the final number density for clusters. 
Carroll, Press and Turner (1992) describe clearly the role 
of the cosmological constant in the rate of structure for- 
mation. As flo is reduced from unity, the rate of growth is 
suppressed, but somewhat less so in the presence of a cosmo- 
logical constant. Thus, in the open case, linear growth stops 
when (1 + z) ~ Qq 1 , when the universe effectively becomes 
curvature dominated, but growth stops more recently, when 

(1 + z) ~ Q 3 , in the flat case, when the universe effectively 
becomes dominated by the cosmological constant. 

There are clear qualitative differences between Einstein- 
de Sitter universes and open models from iV-body simula- 
tions. There is much more dominance of clusters and groups 
of galaxies at earlier times in the open models, which are 
then frozen in; however, open universes do not display so 
prominently the large-scale filaments and other irregular 
structures which occur in Qo — 1 universes. Flat cosmolo- 
gies with a cosmological constant are intermediate between 
these cases. It should be mentioned at this point that the 
results of iV-body simulations require normalisation against 
observations; whilst the foregoing description of structure 
formation rates is not very sensitive to the shape parame- 
ter, r, in the power spectrum, both P and the form of the 
initial conditions (for example, Gaussian or non-Gaussian) 
may effect the resulting abundances on different scales. For 
example, the SCDM model fails to reproduce correctly the 
shape of the galaxy correlation function on scales of tens of 
Megaparsecs, using T = 0.5, but the TCDM model, with 
r = 0.25 does much better. 

We now consider the values of the computed three- 
dimensional shear in each time-slice, and their development 
with time, to see whether comparisons may be made with 
the development of structure in the different cosmologies. 
It is also informative to assess at what point (or in which 
simulation boxes) the maximum contributions to properties 
such as the magnification may occur. 

By simply taking the rms values of specified computed 
components, it is possible, in a simplistic way, to obtain 
values which characterise each time-slice, independent of 
its redshift. This is done (a) before integration along the 
line of sight (which would be necessary to obtain the two- 
dimensional effective lensing potentials, from which the mag- 
nifications are obtained), (b) before conversion to physical 
units (which involve factors to convert from the code units, 
and factors which derive from the evolving box dimensions) , 
and (c) before the application of the appropriate angular 
diameter distance factors (required in the application of 
the multiple lens-plane theory, described in Section 2). Fig- 
ures § |, | and g show the results for the SCDM, TCDM, 
OCDM and LCDM cosmologies respectively. In each, the 
middle curve shows the rms value in each time-slice of the 
sum of the first two diagonal elements of the shear matrix, 
the top curve shows the mean of only the high values of 
these elements, and the lowest curve shows the mean of the 
high values for one of the off-diagonal elements. (The high 
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Figure 4. Curves characterising the time-slices in the SCDM 
cosmology, derived directly from the computed values, and before 
applying any of the conversion factors. The ordinate is in arbitrary 
units. Middle curve: the rms value in each time-slice of the sum 
of the first two diagonal elements of the shear matrix. Top curve: 
the mean of the high values of these elements. Lowest curve: the 
mean of the high values for one of the off-diagonal elements. 
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Figure 5. Curves, as for Figure characterising the time-slices 
in the TCDM cosmology. 

values in each case are those more than la above the mean 
in the rms values.) The growth in the values with time is 
clearly seen in each cosmology. 

In Figures ^| and ^| we have combined the data for the 
different cosmologies. Figure [s] shows the directly computed 
rms values, as above, in the SCDM, OCDM and LCDM cos- 
mologies, and Figure ^ the values in the SCDM, TCDM and 
LCDM cosmologies. (There is some repetition here only to 
make the comparisons clear, and to avoid too much over- 
lapping of the curves.) The plots in Figures |^ and ^ can 
be understood in terms of the discussion at the beginning 
of this section. The Einstein-de Sitter universes, SCDM and 
TCDM, show the most rapid growth of the shear compo- 
nents at late times, reflecting the rapid growth of structure; 
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Figure 6. Curves, as for Figure |4|, characterising the time-slices 
in the OCDM cosmology. 
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Figure 7. Curves, as for Figure characterising the time-slices 
in the LCDM cosmology. 



the OCDM and LCDM results would seem to indicate only 
a limited effect from the cosmological constant. 

Since the real physical dimensions of the simulation 
time-slices evolve with time, it is necessary to introduce fac- 
tors of (1 + z) 3 to the computed values to determine the 
real shear matrix. Doing this appears to dilute considerably 
the effect of structure on the form of the curves, as the ex- 
ample of Figure [To] for the TCDM cosmology shows. (The 
curves are similar in form for all the cosmologies.) The in- 
terpretation of this dilution is that, even though structure 
is forming (to produce greater shear values locally), the real 
expansion of the universe (causing the mean particle separa- 
tion to increase) reduces the shear values in most locations, 
and therefore just outweighs the increases from the forma- 
tion of structure. In this way, the magnitudes of the shear 
component values are seen to reduce slowly with time. 

To obtain the two-dimensional effective lensing poten- 
tials (to which the multiple lens-plane theory may be ap- 
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Figure 8. Curves characterising the time-slices in the SCDM, 
OCDM and LCDM cosmologies, derived directly from the com- 
puted values, and before applying any of the conversion factors. 
The ordinate is in arbitrary units. The curves derive from the 
rms values in each time-slice of the sum of the first two diagonal 
elements of the shear matrix. 



Figure 10. Curves characterising the time-slices in the TCDM 
cosmology, obtained by multiplying the computed values by ( 1 + 
z) 3 in each time-slice. Middle curve: the rms value in each time- 
slice of the sum of the first two diagonal elements of the shear 
matrix. Top curve: the mean of the high values of these elements. 
Lowest curve: the mean of the high values for one of the off- 
diagonal elements. 
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Figure 9. Curves characterising the time-slices in the SCDM, 
TCDM and LCDM cosmologies, derived directly from the com- 
puted values, and before applying any of the conversion factors. 
The ordinate is in arbitrary units. The curves derive from the 
rms values in each time-slice of the sum of the first two diagonal 
elements of the shear matrix. 



plied) it is necessary to integrate the three-dimensional val- 
ues, correctly converted to physical units, along the line of 
sight, and to apply the appropriate angular diameter dis- 
tance factors in accordance with equation This requires 
the integrated three-dimensional computed values (before 
any of the above factors are applied) to be multiplied by the 
factor B(l + z) 2 r d r da /r a , where B = (c/# )(2/c 2 )GM part x 
(comoving box depth) -2 , as described in Section 3.4. When 
the computed values are multiplied by the full conversion 
factors in this way, we see in Figure hll (for the SCDM cos- 
mology, and a source redshift of 3.9) that the peaks are ex- 




Redshift 

Figure 11. The two-dimensional components in each time-slice 
(obtained by integration of the three-dimensional values), con- 
verted to absolute values, including the angular diameter distance 
factors, for sources at z a = 4. Middle curve: the rms value in each 
time-slice of the sum of the diagonal elements of the Jacobian 
matrix. Top curve: the mean of the high values of the summed 
diagonal elements. Lowest curve: the mean of the high values for 
one of the off-diagonal elements of the Jacobian. 



tremely broad, indicating that significant contributions to 
the magnifications and ellipticities can arise in time-slices 
covering a wide range of redshifts, and not just near z = 0.5, 
where the angular diameter distance multiplying factor, R, 
has its peak (for sources at z a — 4). (In this exercise, we 
have used a = 1 for the angular diameter distances. This is 
explained more fully in the next section.) 

The comparisons amongst the cosmologies are interest- 
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Figure 12. The rms values of the sum of the diagonal elements 
of the Jacobian matrix in each time-slice (obtained by integra- 
tion of the three-dimensional values), converted to absolute val- 
ues including the angular diameter distance factors, for sources 
at z 3 = 4. Results are shown for the SCDM, TCDM, OCDM, and 
LCDM cosmologies. 



ing. Figure [l2| shows the integrated rms values (as above) 
multiplied by the full conversion factors and the angular di- 
ameter distance factors for the four cosmologies, assuming a 
source redshift of 4. (The actual source redshifts differ very 
slightly in the different cosmologies.) We see that the LCDM 
cosmology has both the broadest and the highest peak, sug- 
gesting that lenses throughout the broadest redshift range 
are able to contribute significantly to magnifications and 
two-dimensional shearing of images in this cosmology. Be- 
cause of this, the magnitudes of the magnifications are likely 
to be greatest in the LCDM cosmology. Significantly, the 
much higher values for the angular diameter distance factor, 
R, in the LCDM cosmology appear to be the more important 
factor, rather than the existence of structure. The OCDM 
cosmology displays its peak in the rms values at the highest 
redshift, and also has a very broad peak. The TCDM cos- 
mology has both the narrowest range and the lowest peak. 
The differences between the SCDM and TCDM cosmolo- 
gies probably reflect the differences in structure on different 
scales in the two cosmologies, since R is the same for both. 
In general, however, we can say that, for all these cosmolo- 
gies, significant contributions to the magnification and shear 
may arise from lenses at a very wide range of redshifts. 

In a similar study, Premadi, Martel and Matzner 
(1998a) find that the individual contribution due to each 
of their lens-planes is greatest at intermediate redshifts, of 
order z = 1 — 2, for sources located at z s = 5, and Premadi, 
Martel and Matzner (1998b, c) also find very broad peaks 
covering a wide range of intermediate lens-plane redshifts 
for sources at z s =3. 



5.2 Results in the different cosmologies 

In Section 5.1, we assumed a value of a = 1 in the deter- 
mination of the angular diameter distances in the various 
cosmologies. From the output of the shear algorithm one 



is able to obtain an estimate of the dumpiness or smooth- 
ness in each time-slice. Having set the minimum softening 
scale, the code declares the number of particles which are as- 
signed the minimum softening, and one can therefore imme- 
diately obtain the mass fraction contained in clumps, which 
we choose to define by the minimum softening scale. 

For the SCDM cosmology, there is a mass fraction of 
0.026 in clumps in the earliest time-slice at z = 3.6 (next 
to z = 3.9) giving d(z = 3.6) = 0.97, and at z — the 
fraction is 0.17, giving a(z = 0) = 0.83. It is clear that the 
mean value throughout the redshift range is close to 1, and 
almost equivalent to the 'filled beam' approximation. (This 
result is in agreement with Tomita, 1998c, who finds a to 
be close to 1 in all cases.) The multiplying factors, rdTds/rs, 
are very close for the values a = 0.83 and 1.0, although 
somewhat different from the values for a — appropriate 
for an entirely clumpy universe. The fractional discrepancy 
between a — 0.83 and a = 1.0 at the peak of the curves for 
the SCDM and TCDM cosmologies is 5.2% for a source at 
z s = 4, 2.4% for z s — 2, and for sources nearer than z B = 1 
the discrepancy is well below 1%. 

In the TCDM cosmology with shape parameter 0.25, a 
falls to 0.88 at z = 0. The higher value in the TCDM model 
confirms that the SCDM cosmology (with shape parame- 
ter 0.5) has more dumpiness at late times, allowing higher 
values of magnification and shear to occur. The other cos- 
mologies also have high (although somewhat similar) values 
for the smoothness parameter at z = 0. The values are 0.80 
in the OCDM cosmology, and 0.82 in the LCDM cosmol- 
ogy. Table ^ contains the values of a(z — 0), which there- 
fore represent the minimum values for a, and also the ratios 
R(a = 1)/R[a(z = 0)] at the peaks of the curves for the 
different cosmologies. These have all been evaluated for a 
source redshift of 3.6 to enable direct comparisons to be 
made. 

It has been shown by Barber et al. (1999) that the 
weak lensing statistics show only a small sensitivity to the 
smoothness parameter, a, for values between 0.83 and 1 in 
the SCDM cosmology. Furthermore, the minimum values in 
all the cosmologies were always at least 0.8 throughout the 
redshift range, so that the real values are likely to be closer 
to unity in every case. We have therefore chosen to present 
our results on the basis of a = 1 throughout for all the 
cosmologies. This conclusion is validated qualitatively, since 
the variable softening scheme used in the algorithm ensures 
that almost all rays pass entirely through softened mass. 

The source redshifts, z s , we have chosen throughout this 
work are close to 4, 3, 2, 1 and 0.5, and we shall refer to 
the sources in these terms. The actual redshift values var 
slightly for the different cosmologies and appear in Table 

We show in Figure ^| an example of the distributions 
of the magnifications, fj,, in the LCDM cosmology for four 
different source redshifts. For all the source redshifts and all 
the cosmologies, there is a significant range of magnification. 
From these distributions we have computed the values at the 
peaks, ^ip Ca k. Then, since the distributions are asymmetri- 
cal, we have calculated the values, ^i ow and /Hugh, above and 
below which 97^% of all lines of sight fall, and also the rms 
deviations from unity for the magnifications. (These latter 
rms values have been computed only for the lines of sight 
displaying magnifications between /ii ow and /Jhigh because 
three high magnification events in the LCDM cosmology 



© 0000 RAS, MNRAS 000, 000-000 



Gravitational lensing 15 




Magnification, fi 



Figure 13. Probability distributions for the magnification, for 
z s = 3.6, 2.0, 1.0 and 0.5 in the LCDM cosmology. 



z s 




Mpcak 


rms deviation 


Mhigh 


SCDM 










3.9 


0.835 


0.933 


0.115 
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0.101 
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1.0 
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1.181 


0.5 
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TCDM 










3.9 


0.861 


0.959 


0.091 


1.315 


3.0 


0.877 


0.951 


0.081 


1.286 


1.9 


0.904 


0.966 


0.064 


1.228 


1.0 


0.941 


0.972 


0.039 


1.144 
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0.974 


0.986 


0.019 
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4.0 


0.858 


0.919 


0.115 


1.469 


2.9 


0.884 


0.939 


0.115 


1.469 


2.0 


0.915 


0.942 


0.069 


1.283 


1.0 


0.960 


0.972 


0.033 


1.147 


0.5 


0.985 


0.989 


0.013 


1.062 


LCDM 










3.6 


0.789 


0.885 


0.191 


1.850 


2.0 


0.870 


0.934 


0.108 


1.453 


1.0 


0.944 


0.966 


0.045 


1.191 


0.5 


0.981 


0.987 


0.016 


1.070 



Table 4. Various magnification statistics for the different cos- 
mologies, as described in the text. 

would distort the rms values considerably.) All the values 
mentioned are displayed in Table ^. 

It is interesting to compare the distributions in the dif- 
ferent cosmologies. Figure [G|and[l| show the magnification 
distributions for all the cosmologies for source redshifts of 
4 and 1, respectively. The distributions for the magnifica- 
tions (and also the convergence, shear and ellipticities) are 
all broader in the SCDM cosmology when compared with 
the TCDM cosmology, due to its more clumpy character. 
For source redshifts of 4 the OCDM and SCDM cosmolo- 
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Figure 14. The magnification probability distributions for all 
the cosmologies, assuming z s = 4. 




Magnification, [/, 



Figure 15. The magnification probability distributions for all 
the cosmologies, assuming z s = 1. 



gies have very similar distributions even though the angu- 
lar diameter distance multiplying factors are larger in the 
OCDM cosmology. For high source redshifts the magnifica- 
tion distributions are broadest in the LCDM cosmology (and 
the maximum values of the magnification are greatest here), 
but for lower source redshifts the width of the distribution 
is below the SCDM and OCDM cosmologies. 

We plot in Figures and ^ the accumulating number 
of lines of sight having magnifications greater than the ab- 
scissa values. This is done for all the cosmologies for source 
redshifts of 4 and 1 respectively, and clearly shows the dis- 
tinctions at the high magnification end. In particular, the 
LCDM cosmology exhibits a very broad tail for z s — 4, and 
for this reason we show the accumulating magnifications for 
this cosmology in Figure [l8| for source redshifts of 4, 2, 1 
and 0.5. 

In Figure [l9] we show the magnification, fj,, plotted 
against the convergence, k, for z s = 4, again for the LCDM 
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Figure 16. Accumulating number of lines of sight for which the 
magnification is greater than the abscissa value. The plot shows 
the data for all the cosmologies for z s = 4. 



Figure 18. Accumulating number of lines of sight for which the 
magnification is greater than the abscissa value. The plot shows 
the data for the LCDM cosmology for four different source red- 
shifts as indicated. 
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Figure 17. Accumulating number of lines of sight for which the 
magnification is greater than the abscissa value. The plot shows 
the data for all the cosmologies for z 3 = 1. 



Figure 19. LCDM cosmology: /i vs. re for z B = 4 (crosses). The 
continuous line, shown for comparison, represents fi = 1/(1 — re) 2 . 



cosmology. Departures from the curve represented by the 
values of 1/(1 — ft) 2 clearly arise as a result of the presence 
of the term — -y 2 in the denominator of equation and are 
most pronounced at the high ft end, as might be expected. 
This is true for all the source redshifts and all the cosmolo- 
gies. 

We would generally expect the shear, 7, to fluctuate 
strongly for light rays passing through regions of high den- 
sity (high convergence), and we indeed find considerable 
scatter in the shear when plotted against the convergence. 
This would result in different magnification values along 
lines of sight for which the convergence values are the same. 
Figure ^o] (reproduced from Barber et al., 1999) shows the 
result of binning the convergence values in the SCDM cos- 
mology and calculating the average shear in each bin, for 
sources at z s = 4. We see that throughout most of the range 
in ft the average shear increases very slowly, and closely lin- 



early, and we have found similar trends in the other cos- 
mologies. (At the high k end there are too few data points 
to establish accurate average values for 7.) This result sug- 
gests that there may be a trend towards higher mean mag- 
nification values (as ft increases) than would be the case if 
(7) were constant. 

Figures ^l] and show the distributions in the conver- 
gence, ft, primarily responsible for the magnifications for all 
the cosmologies, for z s = 4 and z s = 1 respectively. 

Table |E| shows the rms values for k in the different cos- 
mologies for all the source redshifts, and compares them 
with the rms values for the magnifications. (This time, un- 
like Table ^, the rms values have been computed from all 
lines of sight, rather than just those with magnifications be- 
tween fi\ om and /Xhigh, so we have not included the values for 
source redshifts of z s — 3.6 in the LCDM cosmology where 
three strong lensing events occurred.) 
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Figure 20. SCDM cosmology: Shear vs. convergence for sources 
at z s = 4 (dots), and the average shear (full line) in each of 
the k bins, which shows a slow and nearly linear increase with 
increasing convergence. 
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Figure 21. The probability distributions for the convergence in 
the different cosmologies, assuming z B = 4. 

The distributions in the shear, 7, (defined according to 
equation are broadest, as expected, for the highest source 
redshifts, and, as before, the LCDM cosmology displays the 
broadest distribution for 7 for sources at high redshift. 

The ellipticity, e, in the image of a source is primarily 
produced by the shear, and we show in Figure ^3] the prob- 
ability distributions for e for z 3 = 4 for all the cosmologies. 
Figure shows the accumulating number of lines of sight 
for which the ellipticity is greater than the abscissa value, 
again for z s — 4 and for all the cosmologies. Figures |25| 
and ^ show the same information for z s — 1. 

For z s = 4, the peaks in the ellipticity distributions oc- 
cur at e = 0.075 (SCDM), 0.057 (TCDM), 0.081 (OCDM), 
and 0.111 (LCDM). For z 3 = 1, the corresponding figures 
are e = 0.034 (SCDM), 0.027 (TCDM), 0.021 (OCDM), and 
0.033 (LCDM). Once again, we see that the LCDM cosmol- 
ogy produces the greatest variation and the highest values 



Figure 22. The probability distributions for the convergence in 
the different cosmologies, assuming z s = 1. 
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0.115 


0.047 
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0.073 


0.031 


0.5 
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3.9 


0.126 


0.052 


3.0 


0.111 


0.047 


1.9 


0.088 


0.038 


1.0 


0.056 


0.025 


0.5 


0.027 


0.013 


OCDM 






4.0 


0.245 


0.066 


2.9 


0.176 


0.056 


2.0 


0.123 


0.044 


1.0 


0.060 


0.024 


0.5 


0.026 


0.012 


LCDM 






3.6 


N/A 


N/A 


2.0 


0.314 


0.064 


1.0 


0.120 


0.032 


0.5 


0.030 


0.013 



Table 5. The rms deviations in the magnification and conver- 
gence for the different cosmologies. 

for the highest redshifts, although this is not the case at 
lower redshifts. 

The ellipticity is very closely linear in terms of 7 
throughout most of the range in 7, for all the cosmologies. 
Some scatter occurs because of the factor containing the 
convergence, k, in equation [5b] . 

Finally, the distance-redshift relation, equation |5^, im- 
plies that there may be an effect on the angular diameter 
distances from the shear. Barber et al. (1999) investigated 
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Figure 23. The probability distributions for the ellipticity for all 
the cosmologies, for z s = 4. 



Figure 25. The probability distributions for the ellipticity for all 
the cosmologies, for z s = 1. 
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Figure 24. The accumulating number of lines of sight for which 
the ellipticity is greater than the abscissa value. The plot shows 
the data for all the cosmologies for z s = 4. 



Figure 26. The accumulating number of lines of sight for which 
the ellipticity is greater than the abscissa value. The plot shows 
the data for all the cosmologies for z s = 1. 



this for the SCDM cosmology, and found that whilst the 
maximum effects of shear on the mean magnification values 
may be at least 10%, only 3.2% of the lines of sight were 
affected in this way. Also, they pointed out that the shear 
has an effect in the distance-redshift relation equivalent to 
increasing the effective smoothness parameter, a, and this is 
true for all the cosmologies. By substituting the mean shear 
value (from the SCDM cosmology) determined for sources at 
z a = 0.5 (where the term in a in equation ^ is largest), they 
found the effect on a (and therefore on the angular diame- 
ter distances) to be completely negligible. Furthermore, the 
importance of the effect reduces with redshift, so that our 
decision to ignore the effects of shear in the distance-redshift 
relation is justified. 



6 DETERMINATION OF THE 

COSMOLOGICAL PARAMETERS 

6.1 Weak shear statistics and the density 
parameter 

Jain and Seljak (1997) have given careful consideration to 
the interpretation of observed shear data (from measured 
galaxy ellipticities) by comparing with the analytical results 
of second order perturbation theory. They claim that non- 
linear evolution leads to non-Gaussian effects in the weak 
lensing statistics which are more easily detected in second, 
and higher order, moments. In particular, the probability 
distribution for the three-dimensional density contrast, 

p 



S(r) 



(45) 



is almost indistinguishable in the different cosmologies, 
whereas the probability distribution function for the con- 
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vergence, k, shows different peak amplitudes and different 
dispersions in the different cosmologies. This is because the 
transition to non-linearity in the evolution of structure de- 
pends primarily on the density contrast alone, but the weak 
lensing signal is strongest in those cosmologies which have 
developed structure at optimum redshifts for given source 
positions, and therefore depends more directly on the rate of 
evolution. Specifically, the non-linear evolution of the power 
spectrum introduces non-Gaussianity to the weak lensing 
statistics, whilst the matter density parameter, Qo, (and to 
some degree the vacuum energy density parameter, Ao) de- 
termines the dispersion in the statistics. 

Jain and Seljak (1997) have computed the expected 
skewness, S, in the convergence, 



S = — ((«= - «) }, 

<7 A 



(46) 



where a represents the standard deviation in the distribu- 
tion for ft. As expected, the skewness was greatest on small 
angular scales, and largest for sources at the lowest red- 
shifts. In general, the LCDM cosmology produced the great- 
est skewness, followed by the OCDM, and finally the SCDM 
cosmologies, for sources at low redshift. Even though they 
have assumed an empty beam scenario with a = 0, and cor- 
respondingly different angular diameter distances from our 
work, our own results for the skewness in k are quite con- 
sistent with theirs. We have computed the skewness, not on 
different angular scales, but as a function of redshift for the 
different cosmologies. It is not possible to state a specific 
angular scale for our data, because of the variable softening 
approach in the shear algorithm. Figure ^ confirms that, 
at redshifts less than about 1.5, the LCDM cosmology gives 
rise to the largest skewness, followed by the OCDM cosmol- 
ogy, and finally, the SCDM and TCDM cosmologies. The 
difference between the SCDM and TCDM cosmologies does 
suggest different forms of structure, which result directly 
from the input shape parameter, F, in the power spectrum. 
In broad terms, it is clear that the skewness decreases with 
source redshift, and decreases with the density parameter, 
precisely as expected from perturbation theory. 

Bernardeau, van Waerbeke and Mellier (1997) also use 
perturbation theory (and the empty beam approach) to as- 
sess how the low order moments in the convergence may 
depend on the cosmological parameters. A number of re- 
sults are predicted. They define the moment, S3, (not the 
skewness, defined above) by 



S3 



(* 3 ) 



(47) 



and predict that S3 oc fig for z s ~ 1, or S3 oc Qq for 
z s <C 1, and they predict a slightly weaker dependence on 
Qo for z 3 > 1. We do not find these exact relationships be- 
cause working in the full beam approximation with variable 
smoothing alters the predictions. However, we do find that 
the value of this statistic does decrease with increasing fio- 
Our results for S3 are displayed in Figure Bernardeau, 
van Waerbeke and Mellier (1997) predict significantly less 
dependence on Ao. 

Jain, Seljak and White (1999) have used ray-tracing in 
iV-body simulations, as described in the Introduction, in an 
attempt to evaluate the density parameter from weak lensing 
statistics. Since they make use of the empty beam approxi- 
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Figure 27. The skewness in re as a function of redshift for the dif- 
ferent cosmologies, indicating broadly that the skewness decreases 
with the density parameter, Qq- 
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Figure 28. The statistic S^(k), defined by equation jl?], as a func- 
tion of redshift for the different cosmologies. It is clear that the 
low density cosmologies have higher values of S3 at low redshifts. 



mation, they describe the possibility of many more lines of 
sight passing through voids from z s — 1 in open universes, 
because of the earlier beginning to structure formation. They 
suggest that this means that the probability distribution for 
k in open universes will have its peak close to the minimum 
value. Higher density universes will have the peak well away 
from the minimum, so that the shape of the distribution will 
be less steep. Moreover, they state that this situation only 
applies on small scales, and on such scales this points to the 
minimum value of k being roughly proportional to f2o- A 
better statistic in the empty beam approximation would be 
that ((ft) — K m in) oc Q,q . However, they admit that the result 
'does depend somewhat on the geometry, as the pathlength 
and angular scale differ between open and cosmological con- 
stant models with the same f2o.' More specifically, this result 
will depend crucially on the approximation used, i.e., the 
empty or full beam approximation, and the procedure for 
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Figure 29. re m ; n vs. redshift for the different cosmologies, show- 
ing the order at low redshift expected from the values of Qq. 
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Figure 30. ((re) re m i n ) vs. redshift for the different cosmologies, 
showing the order at low redshift expected from the values of Qq. 



particle smoothing. They find that the measured minimum 
for k in the open cosmology is close to the predicted value, 
but that in the critical, Qo = 1, cosmology the value was far 
from the empty beam value because there were actually no 
completely empty lines of sight. Our full beam approach will 
of course represent the extreme case. It is therefore reassur- 
ing to see that both our «; m i n (Figure ^) and ((k) — K m in) 
(Figure ^) values for the different cosmologies are in the 
order expected from the values of Qo. However, because of 
the entirely different approach from that of Jain, Seljak and 
White (1999), the results are not able to indicate specific 
values for Qq. 

Jain, Seljak and White (1999) also find for the S3 
statistic that S3 (OCDM) > S3 (LCDM) > S 3 (TCDM) > 
S3 (SCDM) at low redshift. This concurs with our findings 
(see Figure ^8|). They quantify the difference between the 
SCDM and TCDM cosmologies in terms of a weak function 



of the shape parameter, which occurs in the expression for 
S3 from perturbation theory. 

They have also explained in some detail how observa- 
tional data may be used to reconstruct the convergence, 
which is only feasible with large fields where the weak shear 
signal is measurable. However, in view of our quite differ- 
ent approaches, which clearly give rise to discrepancies (al- 
though only in terms of degree), care must be taken when 
interpreting observational data in this way. If we had a bet- 
ter understanding of the form, distribution and evolution 
of the dark and luminous matter in the universe, it might 
be possible to produce simulations and weak lensing exper- 
iments in more realistic scenarios. 



6.2 Weak lensing of high-redshift Type la 
Supernovae 

We have seen in Section 5 the significant ranges in magni- 
fications (dependent on the cosmology) which might apply 
to distant sources. In the absence of magnification (or de- 
magnification) from the large scale structure, it would be 
possible to determine the cosmological parameters, Qo and 
Ao, from the departures from linearity in the Hubble dia- 
gram, provided 'standard candle' sources together with good 
calibration were available for measurement at high redshift. 
This is precisely the route taken by a number of authors, 
most importantly Riess et al. (1998) and Perlmutter et al. 
(1999), both of whom have used high-redshift Type la Su- 
pernova data of redshifts up to 0.97. It is evident from our 
results that full account must be taken of the ranges in mag- 
nification for each of the cosmologies, and in particular the 
cosmologies suggested by the high-redshift Type la Super- 
novae results. 

Both groups of workers, i.e., Riess et al. (1998) and 
Perlmutter et al. (1999), point to cosmologies which are close 
to the Qo = 0.3, Ao = 0.7 cosmological simulation we have 
analysed in terms of weak lensing. Consequently, our results 
from this cosmology are of considerable interest for their 
impact on the determination of the cosmological parameters. 
Using the data we report here, Barber (2000) finds that 
true underlying cosmologies having a deceleration parameter 
go = —0.51 + 0.03/ — 0.24 may be interpreted as having go = 
—0.55, from the use of perfect standard candles (without 
intrinsic dispersion), arising purely from the effects of weak 
lensing. This significant dispersion in go (approximately 2a) 
is somewhat larger than that found by Wambsganss et al. 
(1997) based on a cosmology with Qm = 0.4, Qa = 0.6, 
because of our broader magnification distribution at z = 1. 



7 SUMMARY AND CONCLUSIONS 

In the application of the code for the three-dimensional 
shear, we have had to consider what appropriate angular 
diameter distance values should be applied to the data. The 
values have all been calculated numerically from the gener- 
alised beam equations (see section 4) for the different cos- 
mologies. This was done in each case for a smoothness pa- 
rameter 5 = 1. Our variable softening scheme for the parti- 
cles ensures that nearly all rays pass entirely through soft- 
ened mass, and in addition, we found that the minimum 
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value of a was at least 0.8 (at z = 0) in all the cosmolo- 
gies. The differences in the magnification distributions for 
a — 1 and for the minimum value were almost indistin- 
guishable. The effects of shear on the angular diameter dis- 
tances (through changes to the effective value of the smooth- 
ness parameter, a), were found to be completely negligible 
at all redshifts, so we are justified in ignoring them in the 
distance-redshift relation. Furthermore, they are always ad- 
ditive, making the effective value of a even closer to unity. 

We found, in Section 5.f, that the Einstein-de Sitter 
universes showed the most rapid growth in the 'intrinsic' 
shear components at late times, as expected from the growth 
of structure in these cosmologies. This was as a result of 
studying the computed shear values before the application 
of the angular diameter distance factors, and before conver- 
sion to physical units. The OCDM and LCDM cosmologies 
indicated only a limited contribution from the cosmological 
constant in terms of the growth in the shear values, again 
consistent with the expected evolution of structure in these 
cosmologies. 

When the computed shear values were multiplied by 
the full conversion factors appropriate to the integration, 
together with the angular diameter distance factors, the re- 
sulting curves exhibited very broad peaks, indicating that 
significant lensing may result from structure in a wide band 
of redshifts. Significantly, the LCDM cosmology has both 
the broadest and the highest peak, indicating that this cos- 
mology should produce, for example, the broadest range of 
magnifications. This result would appear to come primar- 
ily from the large values of the angular diameter distance 
factors for this cosmology, rather than any considerations 
about the evolution of structure. However, the broader and 
higher peak for the SCDM cosmology, compared with the 
TCDM cosmology, does indicate the differences in structure 
within them, since they both have the same values for the 
angular diameter multiplying factor, R. 

In Section 5.2, we showed the results for the magnifica- 
tion distributions for the different cosmologies for different 
source redshifts, and these are concisely summarised in Ta- 
ble ^. At high redshift, the LCDM cosmology produces the 
highest magnifications, the broadest distribution curves, and 
the lowest peak values. For sources at z s = 3.6 in the LCDM 
cosmology, 97 1% of all lines of sight have magnification val- 
ues up to 1.850. (The maximum magnifications, not quoted 
here, depend on the choice of the minimum softening in the 
code, although the overall distributions are very insensitive 
to the softening.) The rms fluctuations in the magnification 
(about the mean) were as much as 0.191 in this cosmology, 
for sources at z a = 3.6. Even for sources at z 3 — 0.5 there is 
a measurable range of magnifications in all the cosmologies. 

The immediate implication of these results is the likely 
existence of a bias in observed magnitudes of distant objects, 
and a likely dispersion for standard candles (for example, 
Type la Supernova) at high redshift. 

The magnification versus the convergence showed the 
presence of significant shear, and the mean values of 7 in 
small convergence bins pointed to a possible slow linear 
increase in (7) with n. This would result in a trend to- 
wards higher mean magnification values (as k increases) 
than would be the case for constant (7). 

The distributions for the shear, 7, are broadest, in the 
LCDM cosmology for high source redshifts, and there are 



closely linear relationships between the ellipticity and the 
shear in all the cosmologies. This relationship leads again to 
broad distributions in the ellipticity for high-redshift sources 
in the LCDM cosmology. The peak in the ellipticity distri- 
bution for the LCDM cosmology, for z 3 = 3.6, is 0.111, being 
almost twice the value in the TCDM cosmology. 

Jain, Seljak and White (1999) have expressed the pos- 
sibility of determining the value of the density parame- 
ter from the convergence field from weak lensing statistics. 
They show one-point distribution functions for k, assuming 
sources at z s = 1, for all their four cosmologies and us- 
ing different (fixed) smoothing scales in each. They describe 
the increasing non-Gaussianity of the distribution functions 
as the smoothing scale is reduced, and the increasing tail at 
high k. They also describe the shape of the distribution func- 
tions for negative k, which results from the rate of structure 
formation in the different cosmologies, and claim the inter- 
esting conclusion that the minimum value of n is propor- 
tional to the density parameter. We were able to establish 
from our own work that, in broad terms, the skewness in 
decreases with source redshift, and decreases with increas- 
ing density parameter, as expected. Also, the statistic S3 
(defined in equation (47J) was found to decrease with increas- 
ing f2o- However, we were unable to establish more precise 
relationships, because of the use of the full beam approxima- 
tion. The order of the cosmologies for S3 is the same as that 
presented by Jain, Seljak and White (1999). We also found 
that both K m in and ((k) — K m i n ) are in the correct order for 
the different cosmologies in terms of flo, although, again for 
the same reason, it would not be possible to obtain specific 
values for Qo- 

We now, very briefly, make some comparisons with the 
weak lensing results obtained by other authors. This will 
not be exhaustive, because it has been anticipated that our 
results should be different, for a number of reasons. Primary 
amongst these are the following. First, our results were ob- 
tained using the three-dimensional shear code, which allows 
periodicity, the use of the peculiar potential, the net zero 
mean density requirement, and angular diameter distances 
to every evaluation position within each simulation volume. 
Second, many two-dimensional (planar) approaches may suf- 
fer from inadequate convergence to the true limiting values 
for the shear matrix and angular deflections. Third, we have 
introduced a physically realistic variable softening to the 
method, which requires use of the full beam approximation 
for the angular diameter distances, rather than the empty 
cone approximation used by many authors with either point 
masses, small (fixed) softenings, or small pixellation in the 
planes. 

The magnification distributions of Jaroszyriski et al. 
(1990) for the SCDM cosmology do not have mean values 
of 1, and their dispersions in the convergence for sources 
at z s = 1 and z s = 3 are considerably lower than our val- 
ues, with very little evolution with redshift. In our work, 
the peaks of the ellipticity distributions occured at values of 
0.075 (similar to the value for z s = 4) and 0.034 for sources 
at z a = 3 and 1 respectively, and these are somewhat lower 
than the values of 0.095 (z s — 3) and 0.045 (z s = 1) found 
by Jaroszyriski et al. (1990). Rather surprisingly, however, 
their peak values in the distributions for the shear are quite 
similar to ours, especially for sources at z s = 3. 

Wambsganss, Cen and Ostriker (1998) find magnifica- 
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tions up to 100 in the SCDM cosmology, and correspond- 
ingly highly dispersed distributions, very much larger than 
ours for z 3 = 3. The high magnification tail in the distri- 
butions almost certainly derives from the low value of the 
(fixed) softening scale resulting from the 'smearing' of the 
mass distribution in the 10h~~ kpc x 10h~ kpc pixels. 

Similarly, Marri and Ferrara (1998) show wide magnifi- 
cation distributions, and very high maximum values, again 
which occur as a result of using point particles rather than 
smoothed particles. Their procedure (summarised in Sec- 
tion 1.2) is completely different from ours, and, unlike us, 
they find that the SCDM cosmology has the broadest mag- 
nification probability distribution, followed by the LCDM 
cosmology, and finally their HCDM cosmology. In particu- 
lar, we would disagree with their choice of a = 0, which 
is representative of an entirely clumpy universe, as opposed 
to our finding, that the SCDM universe is close to being 
smooth at all epochs. 

Hamana, Martel and Futamase (1999) find that the dis- 
persions in the probability distributions for k, 7 and fi arc 
all greatest for the Einstein-de Sitter cosmology, and are 
very similar for their open and cosmological constant cos- 
mologies. This is in complete contrast to our findings. Their 
magnification and convergence distributions for z s — 1, 2 
and 3 are much broader than ours, although the distribu- 
tions for 7 in the SCDM cosmology are similar. 

The ranges in magnification from Premadi, Martel and 
Matzner (1998a) appear to be rather similar to ours for 
sources at z s — 3 for their three cosmologies, and the widths 
of the distributions are in the same order for the differ- 
ent cosmologies that we find. This is reassuring because, 
although their method relies on two-dimensional projections 
of the simulation boxes, they include many of the essential 
features to which we have drawn attention, for example, an 
assumed periodicity in the matter distribution, randomly 
chosen initial conditions to avoid structure correlations be- 
tween adjacent simulation boxes, the net zero mean den- 
sity requirement, realistic mass profiles for the particles, and 
use of the filled beam approximation with a smoothness pa- 
rameter, a — 1. They show the average shear for a source 
at z s — 5 contributed by each of the lens-planes individu- 
ally, and find that the largest contributions come from those 
planes at intermediate redshift, of order 2=1 - 2. Similarly, 
they find that the lens-planes which contribute most to the 
average magnifications are also located at intermediate red- 
shifts. In terms of the development of structure in the differ- 
ent cosmologies, they find that the lens-planes contributing 
the most shear and magnification are located at larger red- 
shifts for those cosmologies with smaller Qm- The average 
shear for each redshift has been plotted by the authors. They 
find, for all redshifts, that the values in the SCDM cosmol- 
ogy are many times greater than in the LCDM cosmology, 
which is, in turn, greater than in the OCDM cosmology. This 
is not what we find. 

Fluke, Webster and Mortlock (2000), in applying the 
ray bundle method of Fluke, Webster and Mortlock (1999), 
explain clearly the differences between the empty cone and 
full beam approximations. They use the empty cone approx- 
imation, because of their use of an effective fixed physical 
radius for each particle, equal to \/2x the Einstein radius for 
each. This gives rise to magnification probability distribu- 
tions with /i m i n = 1, arising from the use of a — 0, and high 



magnification tails, arising from the small effective radii for 
particles and clusters. They obtain the weak lensing statis- 
tics for the same cosmologies we have used. For z s — 1, the 
distribution in the magnifications for the SCDM cosmology 
is clearly broader than that for the LCDM cosmology, which 
in turn is broader than the OCDM cosmology. This is also 
true of our data, although the order for the cosmologies, in 
our work, is completely altered for higher source redshifts. 
The authors find, as do we, that the most significant differ- 
ences amongst the cosmologies arise directly as a result of 
the optical depth to the source, which is related to the an- 
gular diameter distances. Moreover, there is a suggestion in 
their results (Fluke, 1999) that Kolmogorov-Smirnov tests 
may be unable to distinguish the cosmologies, even random 
distributions of particles, in terms of their magnification dis- 
tributions, if the same set of angular diameter distance fac- 
tors is applied to each. 

This last statement may prove to be extremely impor- 
tant, as it appears from both our own work and that of 
Fluke, Webster and Mortlock (2000), that the angular di- 
ameter distances (or optical depths) are really the deter- 
mining factor for weak lensing statistics. Moreover, our two 
approaches may represent 'limits' for the true lensing be- 
haviour from the mass distribution in the universe. The 
variable softening facility within our algorithm leads natu- 
rally to the assumption that the universe may be described 
in terms of the full beam approximation. This is, however, 
quite different from the assumptions of most other work- 
ers, who frequently use point particles, or a limited form 
of fixed softening, or small pixellation, and therefore use 
the empty beam approximation. The two approaches give 
rise to quite different expectations and results. The most 
obvious differences are the following. First, strong lensing 
can occur with effectively small particles, leading to high 
magnification tails in the probability distributions. Second, 
magnification distributions in the empty cone approxima- 
tion all have /i m in = 1, whilst in the full beam approxima- 
tion, /imin < 1; this may alter the dispersions in the two 
distributions. Third, the mean values for the magnifications 
can be calculated from the respective angular diameter dis- 
tances in the different cosmologies for the empty beam ap- 
proximation; however, the mean values in the full beam ap- 
proximation are always 1. These points make comparisons 
between methods using the different approximations diffi- 
cult. However, it is probable that the universe is neither 
completely smooth nor filled with galactic-mass point-like 
objects. (Subramanian, Cen and Ostriker, 1999, highlight 
our uncertainty in this area, by suggesting that small dense 
masses formed early during hierarchical clustering may per- 
sist to late times, so that real cluster halo structures may 
depend crucially on the detailed dynamics of the dense pock- 
ets.) The resolution of this question together with a better 
understanding of the form, distribution and evolution of the 
dark and luminous matter content of the universe should 
provide a much clearer indication of the likely weak lensing 
statistics in cosmological iV-body simulations of the future. 
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APPENDIX A: GENERALISATION OF THE DYER-ROEDER EQUATION 

In Section 4, we presented the most simple form for the Dyer-Roeder equation (equation [m]) , which can be solved analytically 
for Oo = 1; = 0, and arbitrary a. The general solution for the angular diameter distance between redshifts of Z\ and 22, in 
such a cosmology, is well documented (see, e.g., Schneider et al., 1992) and from this solution the multiplying factors DdDd s /D s 
are easily obtained. However, this solution applies only to cosmologies with zero cosmological constant. We therefore generalised 
the form of the Dyer-Roeder equation to apply to the cosmologies being studied in this work. The following summary of this 
work is likely to be helpful to others working in this field. 

We started from the generalised beam equation, quoted by Linder (1998a and b): 



— + [3 + q(z)](l + z) 1 — + -(l + z ) 2 D 



In this equation, 
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is the deceleration parameter at redshift z, each value of s denotes a content component of the universe (for example, non- 
relativistic matter, radiation, vacuum energy, etc.), and a s and Q 3 (z) represent the smoothness parameter and the density 
parameter respectively, applicable to the component s, and redshift, z. For two-component cosmologies, to which we have 
restricted our work, s 



Al 



reduces to the Dyer- 



for dust and s = — 1 for the vacuum energy. When s — only, equation 
Roeder equation immediately. Also, with two components only, the equation for the deceleration parameter at the present 
day, reduces from equation A2 to the familiar form, go = §O(0) — A(0), where now we have used £1 and A to represent the 



matter and vacuum energy density parameters respectively. Also, the Hubble parameter is, in general, 

x 1/2 



H(z) = H 



o s (o)(i + * 



xi+s) 



(1+^) 2 



For a two-component universe this becomes: 

H{z) = Ho {fi(0)(l + zf + A(0) - [O(0) + A(0) - 1] (1 + z) 2 } 1 ' 2 . 



For two components only, the generalised beam equation (equation Al) is: 
d 2 D . r „ . , NW „ N _\ dD 



dz 2 
in which 



+ [3 + q(z)] (1 



dz + 2 {1 + ^ 



e Da(z)Q,(z) = 0, 



I [Q(0)(l + 2)-2A(0)(l + z)- 2 ] 
q[Z > n(0)z + A(0)(l-M)- 2 -A(0) + l 

0(2) = O(0)(l + zf [H(z)/H ]- 2 , 

H{z) 
H 

and so 



{O(0)(l- 



A(0)-[fi(0) + A(0)-l](l + z) 2 } 



1/2 



0(0) = O(0)(l + zf {O(0)(l + zf + A(0) - [O(0) + A(0) - 1] (1 + zf} 1 . 



To solve equation 
D( Zl , Zl ) = 0, 
and 

dD(z\, z) 



% boundary conditions 



i=(l + *0 



Hjzi) 
Ho 



(A3) 

(A4) 
(A5) 

(A6) 
(A7) 
(A8) 

(A9) 
(A10) 

(All) 



are set, where the second condition is made by considering the form of the Hubble law for a fictitious observer at the redshift 
Z\. 

However, we need values for the angular diameter distances between any arbitrary redshift values (not always based on 
z = 0) in order to construct values for DdDds / D s at all the required evaluation positions. To do this, equation A5 has to be 
generalised further to apply to any arbitrary redshift, and we can do this by changing the variable to 

l + z , 



1 + 21 



(A12) 
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w then corresponds to the re dshi ft of an object as if viewed by an observer at the arbitrary redshift Z\. Then substituting the 
expression for q(z) (equation A6), equation becomes, after some manipulation, 



d 2 D 



dw 



2 (l + *i) 2 



+ < 3 + 



i [Q,(p)x - 2A(0)aT 



I 1 dD 

O(0)(a; - 1) + A(0)x- 2 - A(0)TT j x dw 
O(0)as 



l]x 2 



0. 



(A13) 



(A14) 



fi(0)^ 3 +A(0)-[n(0) + A(0)- 

with boundary conditions, 
L>Oi,zi) = 0, 
and 

djP( f 1,w) \ w =, x = (i + zi)" 1 {n(o)(i + Zl f + A(o) - [n(o) + A(o) - i] (i + Zl fy 1/2 . (Ai5) 

dw k J 

(In equation [A13| we have written x = (1 + w)/(l + z±) for clarity.) 

We have solved this equation numerically for all the cosmologies, checking carefully that the results are the same as the 
analytical values for the Einstien-de Sitter model. To solve it, we made further definitions to simplify the form of the equation. 
First, for clarity, we directly interchanged w and z, and made the following definitions. 

a = (1 + zi) 2 , 

b = fl{0)/(l + Zl ), 

c = 2A(0)(l + 2i) 2 , 

d = 2fi(0)/(l + 

e s2-2fi(0) - 2A(0), 

/=|an(o)/(i + *i), 

g = n(0)/(l + Zl f, 
h = A(0), 

l EE[fi(0)+A(0)-l]/(l + 2l) 2 

and 



j ee (i + Zl )- X {n(o)(i + zi) 3 + a(o) - [n(o) + a(o) - 1] (i + zi) 2 } 1/2 . 

Then the general two-component equation to solve is 



drD 
~dz? 



3 



6(1 + z) - c(l + z) 



d(l + z) + e + c(l + z)- 2 
with boundary conditions 
D{ Zl , Zl ) =0, 
and 

dD(zi,z) _ . 
dz l*=*i- J ' 



a(l + z) 



/(I 



ff (l + z) 3 + ^-i(l + z) : 



-D = 0, 



(A16) 
(A17) 
(A18) 
(A19) 
(A20) 

(A21) 

(A22) 
(A23) 
(A24) 

(A25) 

(A26) 

(A27) 
(A28) 



Figure [l] shows the result of solving equation 
= 3.6, and Figure |2| shows th 
different cosmologies in Table |^. 



A26, with a = 1, in the different cosmologies, for a source redshift of 
= 3.6, and Figure |2| shows the values of rdrd s /r s , also for a = 1. We have tabulated the ratios R(a = 1)/R(a — 0) for the 
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